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It has been proposed that the microwave-induced "zero-resistance" phenomenon, 
observed in a GaAs two-dimensional electron system at low temperatures in moderate 
magnetic fields, results from a state with multiple domains, in which a large local 
electric field E(r) is oriented in different directions. We explore here the questions of 
what may determine the domain arrangement in a given sample, what do the domains 
look like in representative cases, and what may be the consequences of domain-wall 
localization on the macroscopic dc conductance. We consider both effects of sample 
boundaries and effects of disorder, in a simple model, which has a constant Hall 
conductivity, and is characterized by a Lyapunov functional. 

PACS numbers: 

I. INTRODUCTION 

The spontaneous formation of domain patterns is a frequent occurrence in non-equilbirium 
systems, and it has been extensively studied in fluid systems subject to heating or rotation, in 
non-equilibrium crystal growth, in driven interfaces between immiscible fluids, in roughness 
patterns on fracture surfaces, in chemically reacting fluids, in dielectric breakdown, in liquid 
crystal structures, in nonlinear optics, and in motion of granular materials . In some 

cases, domain formation is truly an example of spontaneously broken symmetry; i.e., there 
are two or more domain structures which are equally possible, and selection between them 
is triggered by thermal noise or some other accident of the history of the sample. In other 
cases, small perturbations due to frozen-in disorder, peculiarities at the boundaries, or other 
deviations of the physical system from a symmetric idealized model are responsible. In this 
paper, we discuss a recently discovered system where spontaneous formation of domains is 
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believed to occur and to play an important role in macroscopically observable properties, 
the so-called microwave-induced zero- resistance state in two-dimensional electron systems, 
in the presence of a moderately-strong magnetic field, at low temperatures. 

Experiments on very high mobility two-dimensional electron systems in GaAs heterostruc- 
tures, in fields ranging from about ten to a few hundred millitesla, have shown large changes 
in the dc resistance in the presence of microwave radiation, whose sign depends on the ratio 
Co = uj/uj c between the microwave frequency to and the cyclotron frequency a; c 0|. In par- 
ticular, decreases in the resistance have been observed when Co is somewhat larger than the 
nearest integer. Moreover, in these frequency intervals, if the microwave power is sufficiently 
high, the resistance can drop by several orders of magnitude, perhaps falling below the ex- 
perimental sensitivity, whence the designation Zero Resistance State. We note that in the 
presence of a magnetic field, where there is a nonzero Hall conductance, zero longitudinal 
resistivity is equivalent to zero longitudinal conductivity: the current is perpendicular to 
the electric field. Here, we find it more convenient to emphasize the conductivity rather 
than resistivity. The longitudinal resistance is most directly measured in a Hall bar geom- 
etry, whereas the longitudinal conductance is measured in an annular Corbino geometry, 
illustrated in Fig. |2 below. 

Two distinct microscopic mechanisms for conductivity corrections have been proposed: 
(i) the displacement photocurrent (DP) [5], which is caused by photoexcitation of electrons 
into displaced guiding centers; and (ii) the distribution function (DF) mechanism, which 
involves redistribution of intra Landau level population for large inelastic lifetimesp]. 

Andreev et. al.\$, have noted that irrespective of microscopic details, once the radiation is 
strong enough to render the local conductivity negative, the system as a whole will break into 
domains of photogenerated fields and Hall currents. In the proposed domain phase, motion 
of domain walls can accommodate the external voltage, resulting in a Zero Conductance 
State (ZCS) in the Corbino geometry, or a Zero Resistance State for the Hall bar geometry, 
in apparent agreement with experimental reports^]. 

Nevertheless, many questions remain to be answered. What determines the arrangement 
of domains and domain walls in an actual sample? Is the domain pattern static, or does 
it evolve periodically or chaotically in time? What are the effects of sample boundaries 
and inhomogeneities? What do the domain patterns look like in representative cases? If 
there are favored positions for the domain walls, what effects will this have on the measured 
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electrical conductance? 

We address some of these questions here in the framework of a simple phenomenological 
model. Further details may be found in a previous publication by the current authors. [8] 



II. MODEL. 
A. General form 



In our model, in the presence of the microwave field, there is a local, non-linear relation 
between the dc current j(r) and the local electrostatic field E(r): 

j=j d (E,r)+j"(E,r), (1) 

where j H is the Hall current and j d is the dissipative current, defined by the requirements 

<3 = °%» °% = Wllere 

^(E) ee dji/dEp. (2) 



o^W^dfildEp. (3) 

These equations must be combined with the continuity equation, V • j = —p, where p is the 
charge density, and additional equations that relate the electric field to the charge density. 
In our simplest model, we assume that the Hall current is given by a linear relation 

f = /zxE, (4) 

and the Hall conductance a H is assumed to be a constant, independent of position. As we 
shall see below, there exists a Lyapunov functiona^^, which we can use to determine the 
steady states, the global conductance and stability conditions on domain walls in the strong 
radiation regime. 

The function j d , in general, will depend explicitly on the position r, due, e.g., to inhomo- 
geneities in the 2D electron system, and the direction of ] d may not be perfectly aligned with 
E. Equation ((TJ) is valid above an ultraviolet cutoff, which for the system under consideration 
is probably of the order of the cyclotron radius, of order 1 pm. 
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Writing E = — V0, we may relate changes in the electrostatic potential to changes in 
p through the inverse capacitance matrix W: 

5(j)( r ) = J d 2 rW{r,r')5p{r'). (5) 

If a time-independent steady state is reached, then we have simply V • j = 0, and the precise 
form of W is unimportant. 

In a Corbino geometry, one specifies the potential on the inner and outer boundaries of 
the sample, and one looks for a solution for 0(r) consistent with these boundary conditions. 
If a H is a constant, then the Hall current cannot contribute to V • j in the interior of the 
sample, so it does not appear in Kirchoff's equations. Consequently, the solution for </>(r) is 
independent of a H , and we may, for simplicity, set a H = 0. To recover the Hall current, one 
inserts the solution E into the second term in 

The symmetry condition on a d a ^ allows us to define a Lyapunov functional for j d , as 

r /-e(x) 

G[<f>) - 

A variation of (jOJ is given by 



J d 2 xg(E), 9 = 1* dE'.j d (E>) (6) 



8G= d 2 xV ■ j d 5(f) - / dsn ■ j d 8<P (7) 

J J bound 

The second integral vanishes on equipotential boundaries, or in the absence of external 
currents. The extrema of G are found to be steady states, with V • j = 0. Using the 
positivity of the inverse capacitance matrix W, one may show that G[(j)(t)] is indeed a 
Lyapunov functional, i.e. a non-increasing function of time, so that its minima are stable 
steady states. In general, G may have multiple minima. Any initial choice of 0(r) will 
relax to some local minimum of G, but not necessarily the "ground state" with lowest G. 
Nevertheless, we might expect that the system is most likely to wind up in a state with 
G close to the absolute minimum, and in the presence of noise, we might expect that the 
system might tend to escape from high-lying minima. 

Using the boundary term in (JZJ), the current across a Corbino sample is equal to the 
first derivative of G with respect to the potential difference 5V between two edges, and the 
conductance is given by the stiffness, or the second derivative: 

8I_ J^G_ 

SV d{5V) 2 ' 1 ; 
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B. Strong Radiation. 

We now consider the clean system (E^ = 0) in the regime of strong microwave radiation 
at positive detuning, i.e., with microwave frequencies slightly larger than the cyclotron 
harmonics uj = mu) c , where m — 1, 2, . . .. Both DP and DF mechanisms produce a regime 
of negative longitudinal conductivity around E = 0, which implies a minimum of g{E) 
at a finite field of magnitude E = E c . If the radiation is strong enough, E c can become 
the global minimum, and the Lyapunov density describes a dynamical phase of spontaneous 
photogenerated fields. Then g can be expanded as 

<? (E) = g (E c ) + 1 -{E - E c fa c + A|V • E| 2 , (9) 

where the subscript denotes the absence of disorder. A representative g and the correspond- 
ing dissipative current j are shown in Figs, and [flp. In order to satisfy equipotential 
boundary conditions and the constraint Vx E = 0, domain walls between different directions 
of E must form. The field derivative coefficient A ~ o~ c l^ w implements the ultra-violet cutoff, 
giving a finite domain- wall thickness scale l^ w . Domain walls give a positive contribution 
to G of order a c E^ld w per unit length. In the absence of disorder, the system will simply 
minimize this residual domain wall contribution, subject to aforementioned constraints. 

Local stability requires that g is concave near E c . By rotational symmetry go is a 'Mexican 
hat' function with a flat valley along |E| = E c . Away from domain walls, the fields are 
'marginally' stable. The clean system has zero conductance in the limit when the system 
size is much larger than l^ w . The system can accomodate a potential difference by moving 
domain walls and by changing the direction of E while keeping |E| = E c , and keeping go 
constant everywhere away from domain walls. The current density is maintained at zero 
everywhere, except for possible contributions that vanish in the limit — > 0, which is a 



property of the clean ZCS 



y 



III. DOMAIN STRUCTURE IN THE ABSENCE OF DISORDER 

Although the Lyapunov cost per unit length of a domain wall depends on the precise 
angular difference of the orientations of the electric fields on the two sides of it, the pattern 
of domain walls in the simplest cases, in the absence of disorder, can be understood as an 
attempt to minimize the total domain wall length, subject to the boundary constraints. If 
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the potential is a constant, which we may take equal to zero, everywhere along the boundary, 
we are led to the trial solution: 

0(r) = E c d(r) (10) 

where d(r) is the distance to the nearest boundary. The negative of this solution, with E c 
replaced by —E c , is an equally good solution, with the same Lyaponov cost. Domain walls 
for these solutions occur along curves which are the loci of points that are equidistant from 
two different points on the sample boundary. Domain walls may end, or one domain wall 
may split into two, at special points, which are equidistant from three different points on the 
boundary. Some simple examples, illustrating these possibilities, are given in Figs. El and H3 

Clearly, the trial solution (fTUj) satisfies the requirement that E = E c everywhere except 
on the domain walls. Moreover, the requirement E = E c requires that in the absence of 
singularities, two nearby equipotential contours must be separated from each other by a 
constant distance 5, equal to the difference in the potentials divided by E c . Thus, if one 
starts from a point on the sample boundary, the potential 0(r) must be given by (|10p or its 
negative at least until one encounters a domain wall. If one adds additional domain walls, 
one may construct solutions different from these, but only by increasing the total length 
of the domain walls and presumably increasing the Lyapunov cost. Note that one cannot 
simply displace a domain wall to one side or another without violating the condition E = E c , 
or introducing a discontinuity in (i.e., an infinite electric field), which we do not allow. 

In the circularly symmetric Corbino geometry, illustrated in Fig. El we see that if there is 
no voltage difference between the inner and outer edges, the domain wall must sit at radius 
R = R m where R m is the arithmetic mean of the radii of the inner and outer edges. The 
domain wall has a Lyapunov cost SG = 27tRg w , where g w is the cost per unit length of a 
180-degree domain wall. If the the potential <fi at the outer edge is increased by a small 
amount V, while the potential at the inner edge is held constant, the position R of the 
domain wall must change by an amount 5R = ±V/2E C , where the sign depends on which of 
the two solutions, (fTU|) or its negative, one is starting from. Consequently, since / = dG/dV, 
there should be a small but non-zero photocurrent, / = ±irg w /E c between the outer and 
inner edges, at V = 0. 

Non-zero V breaks the perfect degeneracy between the two steady state solutions, as 
the solution with negative 5R has the smaller value of G, and it has current in the opposite 
direction to the voltage drop. If the system can switch from one state to another, the current 
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will change direction discontinuously when V passes through zero. In the absence of noise, 
or if V is changed sufficiently rapidly, however, the system may be trapped in one state or 
the other, in which case there will be no discontinuity at V = 0. 



IV. LONG WAVELENGTH DISORDER. 



In an inhomogeneous system, there will generally be a non-zero electrostatic field, E^(r) = 
— V0d(r) present in the thermal equilibrium state, with no microwave radiation or bias 
voltage. (Recall that it is the electrochemical potential, which is the sum of the electrostatic 
potential and the internal chemical potential, that is constant in equilibrium). We may ask 
how this disorder field will alter the ZCS. 

At weak disorder field, |E^| << E c , we expand the Lyapunov density in a region where 
E^E C , 

g(E, E d ) = g (E) - a 1 (E) E • E d (r) + 0(E 2 d ), (11) 
which yields a current density 

j"(r) = - gl E, + fl6 "°f - Ed E, (12) 

where X' = dX/dE, and the coefficient o\ depends on the microscopic mechanism. A 
typical g in the presense of disorder is displayed in Fig. For a simple model based on 
the displacement mechanism, we estimate that a c and ai(E c ) should be close to the dark 
conductivity, while a[ <C a c /E c . Other models may have quite different values of a±. 

Our use of (1111) assumes that the length scale of the disorder field is greater than the 



domain wall thickness ld w , which is the same as our short wavelength cutoff. [9] In practice, 
however, we will only be interested in fluctuations on length scales much longer than this. 
In two dimensions, if the disorder potential (pd does not have correlations on a length scale 
much longer than l^ w , and if E^Ui <C E c a c , then the gain in G obtained by aligning E 
with Erf will be too small to overcome the cost of introducing new domain walls, [q] Then 
even in the limit of a very large sample, the number and positions of domain walls will be 
determined by the sample shape and boundary conditions, just as if there were no disorder 
present. By contrast, if the disorder has correlations which fall off slowly at large distances, 
then for the model given by (|TT|) . with a fixed root-mean-square value of E^, it will pay to 
introduce a network of new domain walls, if the sample is sufficiently large, and we would 
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expect that local characteristics of the domain wall pattern will become independent of the 
sample size. 

A. Stability Conditions. 

Stability requires that the local conductivity tensor a d , defined by (J2J), has non-negative 
eigenvalues at every point. The lower (transverse) eigenvalue is given by 

O = ^-Vf-^ + 0{Ed) 2 > ) (13) 

so marginal stability occurs at E = E c + a[E d ■ E/cr c . 

In a steady state, the normal current density across a domain wall must be continuous. 
If n is the normal, and Ei and E2 are the fields on the two sides, we find from (fT3|) and fjllj) 
that 

<r-{Ex) Ex • n = a_(E 2 ) E 2 • n + 0{E 2 d ) (14) 

When the normal component of E has opposite signs on the two sides, (as occurs in the 
clean limit), (jUj) can only be satisfied for a^(Ei) = a^(E 2 ) = O(Ed) 2 . This restricts 
the fields at the domain wall to be at their respective marginally stable values. Among the 
other consequences of these equations, it can be shown that if a path along a set of domain 
walls forms a closed loop, the integral of the "two-dimensional disorder- charge density" 
g d (r) = V ■ E d /27r, over the area enclosed by the loop, must be equal to zero, up to possible 



corrections of order -E^-js] 

To a lowest order in -|^, the field E will have magnitude E c within each domain. Coupled 
with the fact that the parallel component of electric field is continuous across the domain 
wall, this implies that the domain wall between any two domains is oriented in such a way 
as to bisect the angle between the field lines of the two domains. 

Finally, we note that generically, the differential conductance of a sample in the Corbino 
geometry can be obtained by solving for the conductance of a linear system with local 
conductivity given by cr d (E(r)), in series with resistive elements along the domain walls, 
which arise from movement of the domain walls in response to a variation in the applied 
bias V. There could also be discontinuities in the current at discrete values of V, if the 
system jumps discontinuously from one minimum of G to another. We find that for weak 
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long-wavelength disorder, the scale of the macroscopic conductance is set by the domain-wall 
contribution. 



B. One-Dimensional Disorder. 

The simplest example to consider is the case of one-dimensional disorder, where (f) d is 
a function of y, independent of the x-coordinate. At wavelengths much larger than ld W , 
if we do not specify the total voltage drop V , the Lyapunov functional G is minimized if 
the system breaks up into parallel domains, so that E is everywhere aligned with E^, and 
j d = 0. These conditions determine E(y) via (fT2"|) . They are consistent with the boundary 
conditions for a rectangular Corbino geometry, which has periodic boundary conditions at 
x = and x = L, and specified potentials at y — and y = W, provided that the voltage 
difference V satisfies 

rW 

V[I]= / dy{W{y)-El{y)). (15) 
Jo 

Thus, at strong radiation intensity and zero current, we see that that the domain walls form 
precisely at maxima and minima of <ftd, at positions yi, i = 1, ...N, where E d (y) changes sign. 
For a general choice of 4>d{y), when no current is drawn from the sample, there will be a 
non-zero photovoltage, whose value is independent of the magnitude of <fid, to lowest order, 
and is determined by the positions y$ : 

V(j = 0) = E C (^-1) N W - 2 D- 1 )^) + ( 16 ) 

where we have assumed that the odd values of i correspond to maxima in <pd- 

If there is a non-zero current jy (in the y-direction), the domain walls will be displaced, 
leading to a change in voltage proportional to jy. For an appropriate choice of jy, such that 
the voltage resulting from jy cancels the zero-current photovoltage, we can get V = 0. Near 
jy = 0, the differential conductivity a yy is given, to lowest order in cf) d , by 

1 = 2E C 1 _ 2E C 

If P(Ed) is the probability distribution for E^ at a random point, it can be shown that 
E^ 1 = -P(O). If Ed is taken from a Gaussian distribution, then Ed = (2 r K) 1 l 2 E d ms , where 
E d ms is the root-mean square value of E^. If (pd is a single sine wave, then E = 2 1 ^ 2 7iE d ms . 
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The transverse differential conductivity a xx , which can also be calculated using ()13|) and 
(H2J), with j J = 0, is given, to a first order in 0^, by a xx = o~d < \E^\ > /E c , where < \E^\ > 
is the mean value of \Ej\. For a gaussian distribution, one has a xx = (2/ii)a yy , while for a 
single sine wave, a xx = (4/ii 2 )a yy . 

In Fig. Ek, we show potentials <f>(y) for the cases j y = and V = 0, corresponding to a 
particular choice of the one- dimensional disorder potential <pd{y)- 

C. Two-Dimensional Disorder Potentials 

1. Egg-carton potential 
Perhaps the simplest 2D choice for (f> d is the periodic separable egg-carton potential, given 

by 

4>d = Asm(k x x) + B sm{k y y) . (18) 

We consider here a finite system, a rectangle L x by L y , with periodic boudnary conditions 
on the current and the electric fields. This, of course, is equivalent to considering an infinite 
system. For bias voltage V = 0, we construct a trial solution for <p by placing domain walls 
on the lines x = 7r(2m + l)/2k x and y = n(2n + l)/2k y , for integer n and m. Our trial 
solution will have a constant electric field in each rectangular domain: 

E = E c (±x cos 6 ±y sin 6 ) . (19) 

Here, 9$ is the angle between the field and the y axis. Combining equations (fTTJ) and ()19j) . we 
see that, to a linear order, minimization of G is equivalent to maximization of J d 2 xE ■ E^. 
The value of G is minimized with the choice #o = arcta,n(Ak x / Bk y ). The domain walls 
satisfy the neutrality condition that J qd(r)d 2 r = in each domain, and the electric field is 
the gradient of a continuous potential as required. 

In the symmetric case, k x = k y , A = B = E^/(2 1 ^ 2 k x ), the domain walls form a square 
array, and #o = vr/4. We have carried out a numerical minimization of G, and have plotted 
the "two dimensional charge density", q(r) = V-E/27T, for the symmetric egg-carton poten- 
tial, for zero bias voltage in Fig. Efc.. Domain walls, appearing here as line singularities, form 
a square array. Evidently, the numerical solution we obtained is very close to the variatonal 
one above. 
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We now construct a variational solution for the case where there is a nonzero voltage drop 
V along x-direction. Assuming that the electric field is still constant within each domain, 
we allow the centers of domain walls that are perpendicular to the voltage drop to shift 
by amount ±5x and rigidly rotate walls through an angle ±59. The electric fields in the 
various domains will now lie in the directions ± (9q + 89) and ± (tt + 9q — 59), so that the 
domain walls still bisect the angle formed by the field lines in the neighboring domains. In 
Fig. 03d, we plot the (numerically obtained) "two dimensional charge density" , for a bias 
voltage corresponding to an average electric field equal to 0.25E C in the x direction. The 
numerical solution is very much like the variational one just described in that domain walls 
perpendicular to voltage drop are shifted and tilted into a herringbone pattern, while the 
walls along the voltage drop are essentially unchanged. 

The variational solution above can be used to evaluate conductivity a xx . For a voltage 
drop V along ^-direction, we have 

2irV 2n 
— — — = A5x sin 9 ±—59 cos 9 . (20) 

Lyapunov functional G will then be 

G(V) = Go + G x f ± (k x 5x) 2 sin 2 9 + 59 2 ^ + ^ {^j sin 2 9 J J , (21) 
where Gq is the Lyapunov "energy" for the case of V — and Gi = 



i^l^bzhi J (Ak x ) + (Bk y ). Using EqUH we mininimize the quadratic form above. By 
taking derivatives of the resulting expression with respect V, as prescribed by Eq|Hl we 
obtain the following expression for a xx = j^G'^O): 



a, 



Trax ^/{Ak x y + {Bk y f 1 + fi sin2 (jr 



/; l + gsin 2 ^o(^) +fcos^ 



(22) 



For the symmetric case, this reduces to a xx = 0.83aiE®/E c . As in the ID case, the conduc- 
tance is of order <JiE v d ms / E c , reflecting the pinning effect of the "disorder potential" 0^. 

Numerical solutions, referred to above, have been obtained by numerically minimizing 
the Lyapunov functional. We employ a MATLAB-based routine to minimize the functional, 
which we discretize on a triangular lattice (typically containing around 14,000 points) in 
order reduce any perturbations arising from lattice anisotropy. The potential 0(r) is subject 
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periodic boundary conditions in the ^-direction, and is required to change by a chosen bias 
voltage V in the x-direction. Having started with a random initial guess for 0(r), the code 
allows the potential to relax to a local minimum. The robustness of a minimum, i.e. its 
"globality", can be ascertained by starting with a different initial guess and/or giving the 
solution we arrive at random kicks and plugging that in for the next initial guess. As already 
mentioned, for the case of egg-carton potential, we find that the variational solutions above 
are quite close to the ones we arrive at via numerics. 

2. Disordered case 

The numerical routine allows us to obtain solutions for disorder potentials much more 
complicated than the simple egg-carton potential. For example, in Fig. EJo, we display the 
self-consistent potential 0(r) for the case of <pd containing 20 Fourier components, chosen 
from a gaussian distribution with < |0d(k| 2 > independent of k. The disorder potential 4>d 
is also shown. In Fig. Efc, we plot the "two-dimensional charge density", proportional to 
V • E, corresponding to this solution. Domain walls again appear as line singularities in the 
charge. Although one might expect frustration to reduce the conductance in complicated 
potentials, the conductivity in this case was found to be similar to that for an egg-carton 
potential with the same value of E d ms . 

In both the egg carton potential and the more complicated potential of Fig. |3b, there is 
frustration in the minimization of G. The field E is unable to align perfectly with E^(r), 
because of the conflicting requirement of E > E c , which arises from stability, and V x E = 0. 
Because of this, the dissipative current j d is generally non-zero, of order o\E^ leading to 
circulating dissipative currents within each domain. These are indicated by arrows for one 
domain in Fig. Et- The much larger circulating Hall currents, of magnitude a H E c , are not 
shown in the figures. 

V. DISCUSSION. 

In this paper, we have explored the effects of sample boundaries and long- wavelength 
potential disorder, on the domain pattern and conductance of the microwave-induced zero- 
resistance state, within a simple model. In particular, our model, with a constant Hall 
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conductance, has the simplifying feature that it has a Lyapunov functional. In considering 
the effects of disorder, we assumed a linear coupling to the disorder-induced electrostatic field 
Ed. Furthermore, while we consider non-linear effects due to the self-consistent electric field 
E, we do not consider non-linearities due to changes in the local electron- density itself, which 
may result from the electric fields of the domains. (We might expect density fluctuations 
to be small, while potential fluctuations are large, if the system is far from any screening 
conductors.) For more complicated systems, where there is no Lyapunov functional, we 
do not even know whether the system will reach a time-independent steady state, in the 
presence of strong microwave radiation. 

Even for our simple model, with its Lyapunov functional G, there are many questions 
to be answered. What actually is the behavior of a system with two-dimensional disorder, 
in the limit of a large system size? Will there be a large number of metastable steady 
states, as in a glass, with associated hysteresis in the observed current /volt age curve? If we 
assume that for all values of the potential difference V between the inner and outer edge, 
in the Corbino geometry, the system can reach the state with lowest value of G, will the 
conductance be finite in the limit of large system size, with magnitude of disorder held fixed, 
or will it, perhaps, go to zero in this limit? 

Experimental work in the field is also in its infancy. The self-consistent electrostatic 
potential <f>(r) generated by the domains can in principle be measured with minimal distur- 
bance of the sample, using, for example, a single electron transistor device. Such a method 
has been used successfully to explore various aspects of the microscopic properties of the 
integer and fractional quantum Hall effect [lfl| . and we may expect important information 
about domain structures to come from such measurements in the future. 

We remark that Alicea et al.^| have shown that a Lyapunov functional exists, rather 
generally, very close to the dynamical phase transition, where the microwave power is such 
that the longitudinal conductivity first becomes negative at E = 0. However, the form of 
their functional is different from the one employed in the present paper. We are interested 
here in systems under strong microwave irradiation, far from the phase transition, so the 
Lyapunov functional of Alicea et al does not apply. 
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FIG. 1: Model Lyapunov functions g and dissipative current j£ for a uniform electrostatic field E x 
along the x-axis (E y = 0), in the presence of strong microwave irradiation, (a) Lyapunov function 
go(E x ) in the absence of disorder. Operating points are at the minima of g, where E x = ±E C . 
(b) Dissipative current corresponding to go. The dissipative current, parallel to E, vanishes at 
the operating points E x = ±E C . There will be a non-zero Hall current in the y-direction. (c) 
Lyapunov function in the presence of a uniform disorder field in the x-direction. 
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FIG. 2: Domain structure for a circularly-symmetric Corbino geometry, in the macroscopic limit, 
without disorder. Solid circle, midway between the inner and outer edges is the location of the 
domain wall when there is no voltage difference between the two edges. Solid arrows show the 
direction of the electrostatic field E, while dotted arrows show the direction of the circulating Hall 
currents. An alternate solution of the equations has the same position of the domain wall, but 
reversed directions for the electric fields and Hall currents. In a Corbino conductance measurement, 
one applies a voltage difference between the inner and outer edges, and measures the electrical 
current from one edge to the other. 
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FIG. 3: Other possible domain structures in a sample without disorder. Heavy solid lines show- 
domain walls; light solid lines show representative equipotential contours; and double lines show 
sample boundaries. Arrows show direction of electric field E. Panel (a) shows structure for a 
rectangular sample. Panel (b) shows how a domain wall can end inside the sample, in the vicinity 
of a curved corner. Panel (c) shows how a domain wall can split into two. 
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FIG. 4: Examples of a disorder potential 4>d{ v ) and the resulting self-consistent potential 4>{r). 
Different vertical scales were used for <f> and (f>d- Panel (a) shows an example of one-dimensional 
disorder, where <fid and <f> depend only on the coordinate y. Solutions are shown for two different 
boundary conditions, corresponding to zero net voltage drop and to zero longitudinal current flow. 
Panel (b) shows <fi(r) resulting from a particular two-dimensional disorder potential 4>d, indicated 
in the lower plot. The disorder potential contains 20 Fourier components, whose amplitudes were 
chosen from Gaussian distributions with variance independent of the wave-vector k. Periodic 
boundary conditions were imposed on <f> and 4>d- 
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FIG. 5: Gray-scale plots of the "two-dimensional charge density" q(r) = V • E/27T, for three 
examples discussed in the text. Domain walls are visible as line singularities in the charge density. 
Panel (a) shows q(r), in one unit cell, when <fid is the periodic symmetric egg-carton potential, 
(fid cos x + cos 2/) with zero net voltage drop across the sample. The arrows show the direction of 
electric field within each domain. Panel (b) shows q(r) for the same choice of (fid, with a non-zero 
voltage bias, corresponding to an average electric field E = 0.25E C in the x-direction. Once again, 
the arrows show the direction of field within each domain. Panel (c) shows q(r) for the more 
complicated example of <fi and (fid illustrated in Fig. Arrows here show circulating dissipative 
currents in one domain. 



